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Data recorded by the DO experiment at the Fermilab Tevatron Collider are analyzed to search for 
neutral Higgs bosons produced in association with b quarks. This production mode can be enhanced 
in the minimal supersymmetric standard model (MSSM). The search is performed in the three b 
quark channel using multijet triggered events corresponding to an integrated luminosity of 1 fb _1 . 
No statistically significant excess of events with respect to the predicted background is observed and 
limits are set in the MSSM parameter space. 
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Supersymmetry (SUSY) [l| is a popular extension of 
the standard model (SM) which overcomes the hierarchy 
problem associated with electroweak symmetry breaking 
and the Higgs mechanism. In contrast to the SM, where 
only one Higgs doublet is required to break the SU(2) 
symmetry, SUSY requires the presence of at least two 
Higgs doublets. In the MSSM five Higgs bosons remain 
after electroweak symmetry breaking; three neutral: h, 
H, and A - denoted as <f>, and two charged: i? ± . The 
Higgs sector can be parameterized by tan /?, the ratio of 
the two Higgs doublet vacuum expectation values, and 
m,Ai the mass of the pseudo-scalar Higgs boson A. 

The Higgs-quark couplings in the MSSM are propor- 
tional to their SM counterparts, with the exact factor 
depending on the type of quark (up- or down-type) and 
on the type of Higgs boson. For large values of tan/3 
at least two Higgs bosons (either A and h, or A and 
H) have approximately the same mass and couplings 
to down-type quarks, which are enhanced by a factor 
tan j3 relative to the SM ones, while the couplings to 
up-type quarks are suppressed. In this large tan (3 re- 
gion the three Higgs boson couplings follow the sum rule 

.9hbb + Slbb + -9Abb ~ 2 x tan2 P x 5h SM ■ In PP collisions 
at y/s — 1.96 TeV at the Fermilab Tevatron Collider, 
the production of Higgs bosons associated with bottom 
quarks (highest mass down-type quark) is therefore, in 
these cases, enhanced by a factor 2 x tan 2 (3 relative to 
the SM. Due to the tan j3 enhancement, the main decay 
for all these bosons is 4> — > bb (the branching fraction, 
B((f> — > bb), is w 90%). The enhanced production and 
branching ratio make the final state with three b jets an 
important channel in the search for MSSM Higgs bosons 
at large tan/?. At a hadron collider this final state has 
a large background from multijet production which is 
poorly modeled by simulation, making the search for this 
topology very challenging. 

MSSM Higgs boson production has been studied at 
LEP which excluded rrih,A < 93 GeV/c 2 for all tan/3 val- 
ues 0. CDF and DO @,H have extended the MSSM 
Higgs boson searches to higher masses for high tan f3 val- 
ues. The result presented in this Letter supersedes our 
previous published result [f|. In addition to including 
more data, this analysis benefits from improved signal 
and background modeling and an improved limit setting 
procedure, which uses only the shape, and not the nor- 
malization, of the final discriminating variable. 

The DO detector is described in Ref. 0]. Dedicated 
triggers designed to select events with at least three jets 
are used in this analysis. Typical requirements are at 
least two jets with transverse momenta px > 25 GeV/c, 
an additional jet with px > 15 GeV/c, and the pp interac- 
tion vertex is required to be reconstructed well within the 
geometric acceptance of the silicon detector. Algorithms 
for identifying b jets at the trigger level are also employed 
in about 70% of the integrated luminosity used for this 
analysis. After data quality requirements the total data 



sample corresponds to 1.02 ± 0.06 fb~ [8|. 

Signal samples are generated for Higgs boson masses 
from 90-220 GeV/c 2 using the leading order pythia event 
generator [§] to generate associated production of and 
a b quark in the 5-flavor scheme, gb — > <f>b. Weights, cal- 
culated with MCFM [Io| , are applied to the signal samples 
as a function of px and rj of the leading b jet which is not 
from the decay of the Higgs boson, to correct the cross 
section and experimental acceptance to next-to-leading 
order (NLO). Multijet background events from the bb, 
bbj , bbjj, cc, ccj, ccjj, bbcc, and bbbb processes (where j 
denotes a light parton: u, d, s quark or gluon) are gener- 
ated with the ALPGEN ll| event generator. The contri- 
butions from other processes, such as tt, Zbb, and single 
top production, are found to be negligible. The ALPGEN 
samples are processed through pythia for showering and 
hadronization. All samples are then processed through a 
GEANT-based 12j simulation of the DO detector and the 
same reconstruction algorithms as the data. A parame- 
terized trigger simulation is used to model the effects of 
the trigger requirements on the simulated events. 

Jets are reconstructed from energy deposits in 
calorimeter towers using the midpoint cone algorithm (l3l | 
with radius = 0.5. Jet reconstruction and energy scale 
determination are described in detail in Ref. [l4[. All 
calorimeter jets are required to pass a set of quality 
criteria with about 98% efficiency and have at least 
two reconstructed tracks within A7?.(track, jet-axis) = 
V(A?7) 2 + (A^) 2 < 0.5 (where j] is the pseudorapidity 
and ip the azimuthal angle). 

We select signal events by requiring at least three and 
at most five jets with pt > 20 GeV/c and \t]\ < 2.5. 
A neural network (NN) based ^-tagging algorithm fl5j ] . 
with lifetime based information involving the track im- 
pact parameters and secondary vertices as inputs, is used 
to identify b jets. Each event must have at least three jets 
satisfying a tight &-tag NN requirement. The single jet 
6-tagging efficiency is w 50% for a light-jet mistag rate 
of w 0.4%. The events with at least two tight 6-tags are 
also kept and used to model the background. Simulated 
events are weighted based on their tagging and fake rate 
probabilities determined from data. Finally, the trans- 
verse momenta of the two highest px jets which are also 
6-tagged are required to be above 25 GeV/c. To further 
increase the sensitivity, the analysis is split into separate 
three-, four-, and five-jet channels. After the event selec- 
tion 3,224 events remain in the exclusive three-jet sample, 
2,503 and 704 events, respectively, in the four- and five- 
jet sample. The signal efficiencies for Higgs boson masses 
between 100 and 200 GeV/c 2 range from 0.3 - 1.2% in 
the three-jet channel (0.2 - 0.6% and 0.01 - 0.12% in the 
four- and five-jet channels). 

The background composition is determined separately 
for each jet multiplicity. The fractional contribution 04 of 
the ith background process is calculated from equations 
linking the 6-tag efficiency, tj, in an event with the Nj 
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observed events: 



events: 



(1) 



Here, j indicates single, double and triple 6-tagged jets 
in an event for different 6-tag criteria, and iV to t is the 
total number of events. The double ^-tagged sample is 
found to be dominated by bbj while the triple 6-tagged 
sample consists of a mix of s» 50% bbb, as 30% bbj, and 
w 20% bbc+bcc. An alternative method to determine the 
background, based on fitting simulated Ht = ^2 Pt jet 
shape templates to the data, confirms the composition of 
the background. 

For every event the two jet pairs with the largest 
summed transverse momenta are considered as possible 
Higgs boson candidates. To remove discrepancies be- 
tween data and simulation originating from gluon split- 
ting (g — > bb), only jet pairs with A7Z > 1.0 are consid- 
ered in the final analysis. 

The following six variables separate the signal from 
the backgrounds and are well modeled by the simulation: 
the difference in pseudorapidity between the two jets in 
the pair; the azimuthal angular difference between the 
two jets in the pair; the angle between the leading jet 
in the pair and the total momentum of the pair; \p bl — 
P~bi\l\P~hi +P&2I1 the momentum balance in the pair; the 
rapidity of the pair; and the event sphericity. Based on 
these kinematic variables, a likelihood discriminant T>, is 
calculated according to: 



V{x\, ....,x 6 ) 



nti pro*) 



nup? g (xi)+nupr 9 (* 



(2) 



where pf 9 (p^ kg ) refers to the signal (background) 
probability density function (pdf) for variable Xi, and 
(x\, Xq) is the set of measured kinematic variables for 
the jet pair. The pdfs are obtained from triple 6-tagged 
signal and background simulation. Two likelihoods are 
built combining simulated samples in the 90— 130 GeV/c 2 
("Low-mass") and 130 - 220 GeV/c 2 ("High-mass") mass 
ranges, providing discrimination at low and high masses, 
respectively. Studies show that this division of the mass 
range gives the best discrimination. 

Several multijet processes contribute to the back- 
ground and the uncertainty on the cross sections is large. 
The bbb component may also contain a contribution that 
is indistinguishable from a signal and cannot be normal- 
ized from the data. To model the background we there- 
fore rely on a combination of data and simulation. The 
distribution of the expected triple 6-tagged (3Tag) sam- 
ple in the two-dimensional T> and invariant mass (M bb ) 
plane, (T>, Mbb), is obtained from the double b- 

tagged (2Tag) data shape multiplied by the ratio of the 
simulated (MC) shapes of the triple and double tagged 



SlZ s (P,M bb ) 



S^ g (V,M bb ) 
S$*C(V,M bb )' 



^fdata 
>2Tag 



(V,M bb ). (3) 



Many uncertainties affecting the simulation cancel in the 
Figure [T] shows T> for data and back- 



iatl ° S^ s (VM bb ) 

ground for the low-mass likelihood in the three-jet chan- 
nel. 
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FIG. 1: Comparison of the low-mass likelihood distributions 
for the 3Tag data and background (bkg exp.) defined by Eq.|3] 
Every event has two entries. Black crosses refer to data, the 
solid line shows the total background estimate, and the shaded 
region represents the heavy flavor component (bbb, bbc, and 
ccb). The distribution for a Higgs boson of mass 100 GeV/c 2 
is also shown. 



The selection cuts on V, 6-tagging, and number of jet- 
pair combinations per event are optimized by maximizing 
the expected sensitivity. The optimal cuts for the like- 
lihood vary between 0.25 and 0.60 depending on the jet 
multiplicity and Higgs boson mass. The agreement of 
the data and the background expectation is verified in 
a control region where the impact of any Higgs boson 
signal is limited, defined by T> < 0.25. The agreement is 
also verified in the case when no likelihood cut is applied. 
Figure O shows the invariant mass for the optimized high- 
mass likelihood cuts. 

Several sources of systematic uncertainties affecting 

"\(V,Mu,) . 



the background shape through the ratio 3Tag 

Eq. [3] are considered. The dominant uncertainty, due to 
the background composition, is estimated by varying the 
ratio of bbj and bbb events in the sample corresponding 
to the uncertainties from the background composition fit. 
In addition, the smaller uncertainties from the kinematic 
dependence of the 6-tagging, the b jet energy resolution, 
the bbb and bbj kinematics, and finally the trigger-level b- 
tag requirement were also evaluated and included in the 
systematic uncertainty on the background shape. 

The Modified Frequentist method [3] is used to esti- 
mate the compatibility of the data with the background- 
only hypothesis (1 — CL b ) as well as to derive limits at 
the 95% C.L. on the cross section times branching ratio 
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FIG. 2: Invariant mass for the high-mass likelihood region for the exclusive a) three-jet b) four-jet, and c) five-jet channels. 
Black crosses refer to data, the solid line shows the total background estimate, and the shaded region represents the heavy 
flavor component (bbb, bbc, and ccb). The lower panels show the difference between the data and the background expectation. 



as a function of tha- Only the shapes (not the normal- 
ization) of the Mf,b distributions are used to discriminate 
signal from background, assuming the width of 4> to be 
narrow relative to the experimental resolution. Table U 
shows the limits and the 1 — CLb values obtained versus 
the hypothesized Higgs boson mass. The low 1 — CLb 



Mass 


a x B 


a X B 1 


— CL b 


(GeV/c 2 ) 


Exp.(pb) 


Obs.(pb) 


(in %) 


90 


170^2 


184 


39 


100 




128 


38 


110 


71 + 29 
' i -20 


69 


52 


120 


41±« 


34 


73 


130 


28 +1 7 2 


24 


70 


140 


25^ 


22 


60 


160 


17±I 


26 


12 


180 


13t* 


23 


4.4 


200 


9±i 


17 


7.0 


220 


7+3 
' -2 


12 


12 



TABLE I: Cross section limits as a function of Higgs boson 
mass. Columns two and three show the expected (Exp.) and 
observed (Obs.) limits on the cross section times branching 
fraction to bb. The total one-sigma uncertainty on the ex- 
pected limits is also displayed. The last column shows the 
value of 1 — CLi,. 



values around a Higgs mass of 180 GeV/c 2 are due to a 
slight excess over the expected SM background. 

The results of this search can be used to set limits 
on the parameters of the MSSM. As a consequence of 
the enhanced couplings to b quarks at large tan/3 the 
total width of the neutral Higgs bosons also increases 
with tan /3. This can have an impact on our search if the 
width is comparable to or larger than the experimental 
resolution of the reconstructed invariant mass of a di-jet 
system. To take this effect into account, the width of 



the Higgs boson is calculated with feynhiggs [17( and 
included in the simulation as a function of the mass and 
tan /3 by convoluting a relativistic Breit-Wigner function 



with the NLO cross section. In the MSSM the masses 
and couplings of the Higgs bosons depend, in addition 
to tan/3 and itia, on the SUSY parameters through ra- 
diative corrections. Limits on tan/3 as a function of mj 
are derived for two particular scenarios assuming a CP- 
conserving Higgs sector [18J: the m™ ax scenario and the 
no-mixing scenario. Since the results depend consider- 
ably upon the Higgs sector bilinear coupling fi, its two 
possible signs are also probed. 

Figure[3]shows the results obtained in the present anal- 
ysis interpreted in these different MSSM scenarios. Sub- 
stantial areas in the MSSM parameter phase space up to 
masses of 200 GeV/c 2 are excluded. No exclusion can be 
obtained for the m™ ax , [i > scenario. This analysis ex- 
tends the mass range over which the search is performed. 
In addition these results benefit from NN 6-tagging and a 
likelihood discriminant as well as improved modeling and 
a robust limit setting procedure, using only the shape of 
the discriminating variable. 
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FIG. 

m h 



3: 95% C.L. exclusion limits in the (tha, tan 0) plane for 
/i = —200 GeV, and no-mixing, /i — —200 GeV and 
fi = +200 GeV. The exclusions from LEP are also shown Q]. 
The width of cf> is larger than 70% of tua above tan/3 = 100 
H — —200 GeV scenario. 
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